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The present state-of-the-art in cooling mechanical resonators is a version of "sideband" cooling. 
Here we present a method that uses the same configuration as sideband cooling — coupling the 
resonator to be cooled to a second microwave (or optical) auxiliary resonator — but will cool 
significantly colder. This is achieved by varying the strength of the coupling between the two 
resonators over a time on the order of the period of the mechanical resonator. As part of our analysis, 
we also obtain a method for fast, high-fidelity quantum information-transfer between resonators. 
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There is presently a great deal of interest in cooling 
high-frequency micro- and nano-mechanical oscillators to 
their ground states. This interest is due to the need to 
prepare resonators in states with high purity to exploit 
their quantum behavior in future technologies [T[ [5j . The 
key measure of a cooling scheme is the cooling factor, 
which we will denote by / coo i- The cooling factor is the 
ratio of the average number of phonons in the resonator 
at the ambient temperature, n?, to the average number 
of phonons achieved by the cooling method, which we 
will denote by (n) coo i- The present state-of-the-art for 
cooling mechanical resonators is sideband cooling, which 
was originally developed in the context of cooling trapped 
ions [3H5] • This method is a powerful and practical tech- 
nique, able to achieve large cooling factors, and these 
have been demonstrated in the laboratory [6T-|15|. 

In the context of mechanical resonators, sideband cool- 
ing involves coupling the resonator to be cooled (from 
now on the "target" ) to a microwave or optical resonator 
(the "auxiliary" ) whose frequency is sufficiently high that 
it sits in its ground state at the ambient temperature. 
The resonators are coupled together by a linear interac- 
tion, and one that is straightforward to implement ex- 
perimentally. In particular, if we denote the annihilation 
operators for the target and auxiliary resonator by a and 
6, respectively, then the full Hamiltonian of the two res- 
onators is 

H = Huiar a + hflb^b + g cos(vt)x a Xb, (1) 

where x a = a + a) and Xf, = b + w are the position 
operators of the respective resonators. The coupling is 
modulated at the difference frequency between the res- 
onators, v = Q — lo. This converts the high frequency 
of the auxiliary resonator so that the two resonators are 
effectively on-resonance, and thus exchange energy at the 
coupling rate g. With this frequency conversion, the aux- 
iliary constitutes a source of essentially zero entropy (and 
thus zero temperature) for the target resonator |16| . 



When the rate of the coupling, g, is significantly 
smaller than the frequency u) of the target resonator 
(so that one is within the rotating-wave approximation 
(RWA) — see, e.g. [T7]), then the linear coupling between 
the resonators is merely excitation (phonon/photon) ex- 
change between the two. If the auxiliary is now damped 
sufficiently rapidly, then the excitation exchange, com- 
bined with the relatively fast damping of the auxiliary 
at effectively zero temperature, extracts the phonons out 
of the target. The cooling factor is merely the ratio of 
the phonon extraction rate to that of the heating rate, 
the latter being the product of the damping rate of the 
target, 7, multiplied by the average number of phonons 
it has at the ambient temperature, ny. It is not always 
possible to determine the extraction rate analytically, but 
if we denote it by r coo i, then 

(n) C ooi = j^- = ( tt~~ ) nr. (2) 

/cool \ ' cool / 

Note that r coo i depends on how fast the auxiliary can 
extract energy from the target, or the ratio g/n. For 
sideband cooling, the RWA requires g/n <C cj/k, and 
thus limits the cooling factor. 

Here we demonstrate that one can cool significantly 
better than traditional sideband cooling by using quan- 
tum control to go beyond the RWA, into the ultra-strong 
coupling regime g ~ w. We first show that a particular 
time-dependence of the coupling rate, g(t), can achieve 
a high-fidelity transfer of quantum states between the 
target and auxiliary resonators within a single resonator 
period. As pointed out in [18], "state-swapping" is one 
way to achieve cooling, as this process will load the 
cold state of the auxiliary into the target. In fact, the 
phonon/photon exchange of the RWA implements state- 
swapping in a time of n/(2g) |19) . However it was shown 
in [18 that using this to cool (which means running tradi- 
tional sideband cooling, but now only for a single swap- 
time) is little better than the usual approach. In con- 
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trast, we show here that numerically optimized control 
sequences will achieve significantly better cooling factors. 
Because this technique requires a relatively small modifi- 
cation to the existing sideband cooling scheme, it can be 
readily implemented with present technology. Further, 
our method not only provides better cooling, but makes 
this achievable over a much wider range of values of the 
auxiliary damping rate, k. While our method always per- 
forms at least as well as sideband cooling for any value 
of g, to obtain the lowest achievable teperatures one does 
require ultra-strong coupling [g ~ oj). This is neverthe- 
less timely, because very recent experiments have demon- 
strated coupling not far from this regime [13J . We note 
that Machnes et al. [20] have previously devised a way to 
go beyond the RWA in the context of trapped-ion experi- 
ments, in which the auxiliary system is a qubit. However, 
their method is not feasible for nano-resonators, certainly 
with present technology, because it requires g ^> oj |13j . 

To begin our analysis we first consider the problem of 
engineering a fast, high-fidelity state-swap between two 
linearly coupled resonators, as this is an important prob- 
lem in its own right. Fast operations on quantum in- 
formation are important due to the the ever present ef- 
fects of decoherence. To obtain such a state-swap, and 
thus an efficient energy transfer without the RWA, we 
examine the algebra generated by the linear coupling in 
conjunction with the free Hamiltonians of the resonators. 
The algebra of these three Hamiltonians suggests that it 
should be possible to engineer a perfect state-transfer op- 
eration between the two resonators, by concatenating the 
evolutions generated by the Hamiltonians in a process of 
"quantum control" [3T]. Up to local operations on each 
resonator, such a concatenation is equivalent to varying 
the coupling g with time. This would allow us to obtain 
efficient energy transfer when g ~ oj, not only achieving 
faster state-swapping, but also better cooling. 

To explore the above conjecture, we simulate the evo- 
lution given by the Hamiltonian in Eq. , in which g is a 
function of time. Since is typically much greater than 
oj (by a factor of at least 100), it is a good approximation 
to assume that the frequency conversion is exact, and set 
f2 = oj and v = 0. The corrections to this approximation 
are of the order of (oj/tt) 2 . (This is, in fact, an RWA for 
the frequency £1, which is distinct from the RWA for the 
target frequency oj, required by sideband cooling.) We 
prepare the target resonator in a state that is confined 
to the space spanned by the 12 lowest Fock states, but 
otherwise completely mixed on that space. The auxil- 
iary is prepared in the ground state, and the resonators 
evolved for a specified time. This allows us to determine 
the quality of the swap merely by calculating the pu- 
rity of the final density matrix for the target resonator. 
If this state is pure, then the evolution has successfully 
transferred all the quantum information to the auxiliary 
resonator. We evolve for a single period of the target res- 
onator, and dividing this time into five equal intervals of 



duration At, we parametrize g(t) by making it piecewise- 
constant on these intervals. Finally we perform a nu- 
merical optimization, using a Quasi-Newton line search 
method |22| , to determine the five piecewise-constant val- 
ues for g(t). For the simulation we use the basis of Fock 
states, including the lowest 25 states for each resonator. 
This achieves an essentially perfect state-swap (a final 
purity of 0.999977) with the following five values of g/oj: 
(1.78,1.45,2.44,1.61,0.195). As a second example, we 
find that a state-swap with a purity of 0.999991 can be 
obtained in 0.7 of the resonators period, with the values 
(2.76,0.474,3.73,0.78,2.59). 

The above results show that, in the absence of deco- 
herence, state-swapping in less than one period is within 
the "control space" of the linear coupling. But this does 
not tell us how well we can transfer the cold auxiliary 
state to the hot resonator in the presence of damping. 
Damping interferes with the swapping process via the 
quantum Zeno. Damping is equivalent to a continuous 
measurement process [HI [23] , and this inhibits the trans- 
fer of energy to the auxiliary. We must therefore simulate 
the optimized cooling in the presence of damping, but 
it is impractical to do this with the simulation method 
used above, as the size of the required superoperators 
is too large. Fortunately in the case of cooling we are 
only interested in the average phonon number, given by 
(n) = (a' a), which is a second moment of operators a 
and . Because the dynamics of the resonators is linear 
(that is, the evolution can be described by a set of lin- 
ear quantum Langevin equations 24-26 ) one can derive 
a closed set of equations for the variances and covari- 
ances of the annihilation operators. Because the means 
of these operators are zero in thermal states, and remain 
zero during the evolution, the covariances are equal to 
the second moments. 

If we define the vector x = (a, (r , b, b^Y, then the ma- 



trix of covariances is C : 
of motion for C is 



(xx') — (x)(x t ). The equation 



C = AC + CA t + G, 



where 
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and n aux is the initial average number of phonons in the 
auxiliary resonator. 

We now use the above equations to perform a numer- 
ical optimization as before, but this time to determine 
the function g(t) that gives the minimum value of (a^a) 
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FIG. 1. (Color online) The average phonon number, (n) coo i, achieved by the cooling method presented here, compared to 
sideband cooling, as a function of the damping rate of the auxiliary resonator, k, and for four values of -jut- The values 
obtained by sideband cooling, for the optimal choice of coupling g, are given by the solid circles/dashed line (red online). The 
lower curves, consisting of various symbols, give the values achieved by the present cooling method, in which g is varied with 
time. The various symbols indicate different values of the total time, r, over which the cooling is achieved. In the following, 
r is in units of oscillator period, (a) From left to right: Star: r = 0.5; Diamonds: r = 0.6; Squares: t = 1; Circle: r = 1.5; 
Diamond: r = 2; Star: r = 3; Squares: r = 5.3; Circles: r = 6. (b) Diamonds: r = 0.7; Squares: r = 1; Circles: r = 6; (c) 
Stars: r = 0.6; Diamonds: r = 0.8; Squares: r = 1; Circles: r = 2; (d) Diamonds: r = 1; Squares: r = 2; Circles: r = 6. 



after a fixed time interval. To do this we take the same 
approach as above, dividing the total evolution time into 
N equal intervals, and making g piecewise constant on 
those intervals. We wish to determine the optimal cool- 
ing over a broad range of the relevant parameters, and 
compare this to sideband cooling. The important param- 
eters are the damping rates of the target and auxiliary 
(respectively 7 and k) , and the average number of ther- 
mal photons in the target at temperature T, tit- We 
note that the average number of photons in the auxiliary 
the ambient temperature, n aux , can be made very small 
with present technology. For example, a 10 GHz stripline 
resonator at 50 mK has n aux = 6.7 x 10~ 5 . The same res- 
onator at 100 mK has n aux = 8.3 x 10~ 3 . We expect n aux 
only to be significant if we are able to cool the target so 
that (a' a) is close to rt aux , and we verify this below. 

By examining the equations of motion, we see that so 
long as riT 3> 1, the evolution, and thus the cooling, 
depends only on the product of 7 and n?, rather than 
each separately. Since S> 1 is the relevant regime for 
present experiments, we therefore need to determine the 
optimal cooling as a function of k and "/nx- Finally, it 
turns out that there is not only an optimal g(t) for each 
value of k/lj, but also an optimal time over which to 
perform the control. We could optimize g(t) and the total 
time in a single shot, but we found that this extended the 
optimization time considerably. We therefore perform 
the optimization for a set of fixed total times, and plot 
the best results for each value of k. 

We now perform the optimization over g(t), with 
ri aux = 0, and plot the results in Figure [T] It is con- 
venient to refer to 7 and k in units of u>, and thus we 



plot the achieved value of (a) a) as a function of k/lj, 
and for four values of the product (7/0;) n^. In Figure [I] 
we also plot, for direct comparison, the values of (a^a) 
that are achieved using sideband cooling. In agreement 
with [5], we find that sideband cooling achieves its best 
performance when n is in the range O.lw - uj. 

We see from Fig. [I] that our "optimal control" cool- 
ing scheme is superior to sideband cooling when n is 
less than the value for which sideband cooling achieves 
its best performance; above that value of k their per- 
formance is very similar. The second key result is that 
the improvement provided by optimal control steadily in- 
creases as the product 777,7^ decreases with respect to u>. 
For (j/oS)nT = 10~ 4 , 10~ 3 , and 10~ 2 , the smallest values 
we obtained for (a) a) are better than sideband cooling 
by factors of approximately 12, 6, and 3, respectively 
(see Figs. [ija)-(c)). For ( r y/uj)nT — 1, while our method 
no longer achieves a lower temperature than sideband 
cooling, it achieves this temperature for a wide range of 
values of k, whereas sideband cooling does so only for a 
single value. 

Almost all the cooling results in Fig. [T] are ob- 
tained with no more than 20 time-segments (that is, 20 
piecewise-constant values for g(t)) per period. The ex- 
ception is the diamonds in Fig. |T|(b) , in which we used 30 
segments in 0.7 of a period. In many case we find that 
between 5 and 10 time-segments per period is sufficient 
to obtain optimal cooling. In Fig. [2] we show the explicit 
control for g(t), for a chosen value of k in Fig. [T^,, along 
with the evolution of (a'a) for this sequence. While we 
have chosen to parametrize g(t) using a piecewise con- 
stant function, one can just as easily use a Fourier series, 
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FIG. 2. (Color online) Here we show the optimal 10-segment 
pulse shape for the interaction rate g(t) (inset), along with the 
time-dependence of the average number of photons, for one 
of the points in Fig. [TJa. In this case n = 100, 7 = le — Qui, 
K = 1.35 x 10 _3 cj, and the duration of of the control pulse is 
0.6 periods of the resonator. 



which removes the need for sharp switching between the 
segments. When optimality is achieved with 20 segments 
per period, this indicates that the highest frequency com- 
ponent required for g in such a Fourier series is no more 
than 20w. For a 50 MHz resonator, this is l GHz, and 
thus experimentally quite practical. We note that as k 
increases, the time, r, required for optimal cooling also 
increases. For k = cj/1000, r ~ 0.7 of the resonator 
period, whereas when k = u)/2, r is typically about six 
periods (the specific cooling times are given in the cap- 
tion of Fig.[l]). 

We now determine the effect of thermal photons in the 
auxiliary resonator (n aux > 0). We perform the opti- 
mization again for a case in which we obtain the lowest 
temperatures: 7 = 10 _6 cl>, tit = 100, and five values of k 
equally spaced on a log scale from 10 -4 w to 10~ 3 u;. First 
we set n aux = and obtain the following set of cooled 
values for (a^a): 10~ 4 x (2.9,4.0,4.3,5.3,7.0). Now per- 



10" 



the new set 



forming the optimization with n aux 
of cooled values for (a^a) arc 10~ 4 x (4.1, 4.4, 5.0, 6.1, 8.4). 
The increase in the average phonon number is approxi- 
mately the addition of n aux . This confirms our intuition 
that the effect of thermal photons in the auxiliary is only 
significant when the target is cooled close to n aux . 

Finally, we consider coupling two auxiliary oscillators 
to the target resonator, in which both coupling rates are 
independent functions of time. Performing an optimiza- 
tion in this case, we find that this configuration does 
not provide any significant improvement in the achiev- 
able cooling. This suggests that the cooling method we 
have presented gives the best cooling that is possible us- 
ing linear elements. 
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